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The reaction ir~p — > lutt~p, lu — » tt + tt~tv° has been studied at 18 GeV/c. The luit~ mass spectrum 
is found to be dominated by the 6i(1235). Partial Wave Analysis (PWA) shows that 6i production is 
dominated by natural parity exchange. The S-wave and 15- wave amplitudes for 6i(1235) — > lutv have 
been determined, and it is found that the amplitude ratio, \D/S\ = 0.269 ± (0.009) s tat ± (0.01) sys 
and the phase difference, <j>(D - S) = 10.54° ± (2.4°) stat ± (3.9) sys . 



I. INTRODUCTION 



According to QCD mesons are bound states of quarks, 
anti-quarks, and gluons. The interaction between quarks 
and gluons is most conveniently simulated by a two part 
potential. The short range behavior of the interaction is 
dominated by the one gluon exchange 'Coulombic' poten- 
tial, while the long range part is dominated by a 'linear' 
confinement potential which is attributed to a collective 
multi-gluon exchange, often modeled as a flux-tube. Me- 
son decays, for example, A— >B+C, necessarily arise from 
the confinement interaction and the corresponding break- 
ing of the flux-tube with the creation of a qq pair in a 
certain 2s+1 Lj state. Several different calculations ex- 



ist in the literature 
pair creation || [I| 



2[ | assumine 
Pq pair creation g, J 



i S 1 



nozarm@jlab.org; Present address: Thomas Jefferson National 



Accelerator Facility, Newport News, Virginia 23606 

» Present address: Department of Physics, Florida State University, 

Tallahassee, FL 32306 

* Present address: Department of Physics, University of Maryland, 

College Park, MD 20742 

§ Permanent address: Rafael, Haifa, Israel 

^Present address: TRIUMF, Vancouver, B.C., V6T 2A3, Canada 
** Present address: Department of Physics, Idaho State University, 
Pocatello, ID 83209 

ttpresent address: Laboratory for Nuclear Studies, Cornell Uni- 
versity, Ithaca, NY 14853 

^Permanent address: Institute of High Energy Physics, Bejing, 
China 

§§Deceased, Dec. 2001 



and with and without final state interactions M Iq, 7\. 
It was pointed out by Kokoski and Isgur j|] that accu- 
rate measurements of amplitude ratios in specific decays 
such as the D-wave to S-wave amplitude ratio, (D/S), 
in b\ — > lutt , cl\ — ► p7r, and h\ — » pw, or the P-wave to 
F-wave amplitude ratio, (P/F), in Ti2 — > fm, could pro- 
vide sensitive measures which distinguish between differ- 
ent model predictions. This is so because these decay 
amplitude ratios are more sensitive to the decay dynam- 
ics than to the hadronic structures (5) . 

The experimental measurements compiled by the Par- 
ticle Data Group (PDG2000) |), yield a wide range of 
\D/S\ values for 6i(1235) — > lott and there are no pre- 
vious measurements of <f>{D — S). The PDG2000 recom- 
mended \D/S\ value of 0.29±0.04 is taken as the weighted 
average of several measurements. The two most recent 
measurements, 0.23 ± 0.03 and 0.45 ± 0.04, reported by 
the Crystal Barrel collaboration from the analyses of the 
uiriTT p] and W7r°7r° ]To[ | final states in pp annihilation at 
rest, are not consistent with each other. The remaining 
experimental measurements from b\ production in 7p pd| ] 
and 7rp jl2[ |l3|, [fj| reactions were based on very low 
statistics and none were sensitive enough to determine 
cb(D-S). 

We have made a high statistics measurement of 
&i(1235) — » lott produced in the reaction ir~p — ► toix^p 
at 18 GeV/c, with u> -> tt+tt^tt and tt° -> 77. The 
final state particles 7r + ,27r _ , and 27 were detected; the 
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recoil proton trajectory was also measured. The uj — > 3n 
decay matrix element was used to estimate the signal 
from the background of other 37r production processes. 
The J pc = l" 1 6i(1235) signal was identified by a Partial 
Wave Analysis (PWA) of the data. The b\ — > (ujtt)s and 
61 — > (ujtt)d amplitudes were determined, yielding new, 
accurate measurements of both \D/S\ and 4>(D — S). 



II. DATA SELECTION AND MAIN FEATURES 
OF THE DATA 

The present measurements were part of the 
Brookhaven National Lab experiment E852, per- 
formed at the Multi-Particle Spectrometer (MPS) 
facility which has been described in detail elsewhere [p| . 
An 18.3 GeV/c ir~ beam, delivered by the Alternating 
Gradient Synchrotron, was incident on a 30 cm liquid 
Hydrogen target. The MPS was equipped with a 4-layer 
cylindrical wire chamber (TCYL) Q for trig germs 
and detection of charged recoil particles, a 198-block 
cylindrical Thallium doped Cesium Iodide detector Q 
around TCYL to veto soft photons, a 3045-element 
lead glass calorimeter (LGD) jl9) for detection of 
photons, and a downstream 2-plane drift chamber, 
located directly in front of LGD, for tagging charged 
particles incident on the LGD. Three proportional wire 
chambers (PWCs) were interspersed between six 7-plane 
drift chambers [EQ| inside the MPS magnet. The first 
two PWCs provided a forward multiplicity trigger. A 
multilayer lead/scintillator sandwich in the shape of a 
picture frame consisting of four box shaped sections 
was placed after the second PWC. This detector in 
conjunction with four scintillators allowed for rejection 
of downstream wide-angle photons that fell outside the 
acceptance of the LGD. 

The trigger for the reaction n~p — > ujtt~p, uj — > 
7T + 7r~7r°, 7T° — > 77 required three forward-going charged 
particles, as well as one large angle charged recoil (i.e. 
the final state proton) in TCYL. A total of 265 million 
such triggers were recorded during the 1995 running pe- 
riod of E852. Photons from ir° — ► 77 were detected in 
the LGD. After requiring charge and energy-momentum 
conservation in addition to topological and fiducial cuts, 
8.2 million events of the type tt + ir~ n~ ^/j and a missing 
mass around the mass of the proton remained. A 2 con- 
straint kinematical fit, requiring a proton recoil at the 
main vertex and a it from the 27's, with a confidence 
level (c.l.)> 5%, resulted in 1.2 million ir + ir~Tr Q Tc~p 
events. An additional kinematical constraint requiring 
the 7T + 7r _ 7r° mass to be consistent with the uj mass was 
imposed. Events with a c.l. > 5% were selected, resulting 
in a final sample of 224 thousand ujir~p exclusive events 
which were then subjected to more detailed analysis. 

Figure [i] shows the 7r + 7r~7r° invariant mass spectrum. 
There are two entries per event, corresponding to the 
two neutral three pion combinations. The hatched re- 
gion corresponds to a cut around the uj peak, defined as 



0.760 < m(7r + 7T"7r°) < 0.845 GeV/c 2 , which dominates 
the spectrum. There is also a significant number of events 
in the region of the oi(1260) and the <22(1320). 




m(7T + 7T 7T°) GeV 

FIG. 1: Invariant mass spectrum for 7r + 7r~7r° 

The 7r + 7r~7r°7r~ effective mass spectrum is shown in 
Figure ||, before (un-hatched) and after (hatched) the 
uj selection. The An spectrum shows two distinct peaks 
after the uj cut, one around the &i(1235) mass, and one 
at the /9 3 (1690) mass. As seen in Figure 0, there is a 
significant, approximately linearly increasing background 
under the uj peak. Its magnitude is on the order of 25%, 
and it arises mainly from the a\/a2 resonances, which 
also decay to 37r. This is confirmed from the PWA results 
discussed in section 6. 

The distribution of four-momentum transfer squared, 
— t, is shown in Figure ||. For the PWA only data with 
0.1 < — t < 1.5 (GeV) 2 were used. In this region, the —t 
distribution was fitted to a function of the form f(t) = 
Pi e P2 * + P 3 e Pi *, with the coefficients Pi = 9.9 ± 0.01 
and P 3 = 7.7±0.10, and the slope values P 2 = 4.5±0.052 
(GeV 2 )" 1 and P 4 = 1.7 ±0.070 (GeV 2 )- 1 . 



III. PWA OF THE lutt- SYSTEM 

The details of the E852 PWA formalism are discussed 
in 1 21, ^2j and a general description on the implementa- 
tion of the formalism, using a PWA code, is given in [E3|. 



In the formalism, the interaction process is divided into 
two parts, the production and the decay. Calculation of 
the decay amplitudes in our PWA formalism follows the 
framework of the isobar model [ p4| , where the decays at 
each vertex proceed through two-body modes. An excep- 
tion to this assumption is the treatment of the uj decay, 
where a direct three-body decay is used. 

In the PWA, each possible resonance is characterized 
by a partial wave, labeled by the quantum numbers 
J pc [isobar\Lrrf , where J pc are the total spin, parity and 
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1.5 2 2.5 3 

m(7T + TT~7T°7T~) GeV 

FIG. 2: Invariant mass spectrum for 7r + 7r~7r°7r~ and uin~ 
combinations. The hatched regions correspond to an u) mass 
cut. 



C-parity of the partial wave, L is the orbital angular mo- 
mentum between the decay products, m is the absolute 
value of the spin projection of the resonance along the 
quantization axis (chosen to be in the beam direction), 
and e is the reflectivity of the partial wave. The reflection 
operator is defined as a rotation about the y-axis by tt 
radians followed by the parity operation. In our analysis, 
the y-axis is chosen in the direction of the normal to the 
production plane. 

The amplitudes, expressed as eigenstates of the reflec- 
tion operator, are constructed from the helicity states, 
to account for parity conservation in the production [^5) . 
There are two advantages to using the reflectivity basis. 
Firstly, the states of different e do not interfere with each 
other, and secondly, in 7rp reactions, there is a direct cor- 
relation between e and the naturality of the exchanged 
particle. 

The total intensity distribution is written as the sum 
of 2k intensities where each intensity is the square of the 
sum of N e interfering amplitudes, and the factor of 2 is 
for e = +1 and e = — 1. k corresponds to the possi- 
ble configurations of the spin at the baryon vertex, i.e. 
spin flip and spin non-flip {k is therefore the rank of the 
spin-density matrix). N e is the number of partial waves 
included in the fit, in a given reflectivity. The intensity 
distribution in terms of production amplitudes, V, and 
decay amplitudes, A, is given by: 
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FIG. 3: The negative of the four-momentum transfer squared, 
—t distribution for the u>tv~ p final state. The distribution 
is fitted to the sum of two exponential functions. The fit 
results yield two slopes, P2 = 4.5±0.052 (GeV 2 )" 1 and P4 = 
1.7 ±0.070 (GeV 2 )" 1 . The arrows indicate the region chosen 
in the PWA: (0.1 < -t < 1.5 (GeV) 2 ). 
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The subscript a denotes a set of parameters specifying 
the interfering amplitudes, for instance, total spin of the 
state, </, its parity, p, the component of the total spin 
along z, m, and the orbital angular momentum between 
its decay products, L. t is a set of independent variables 
which specify the configuration of the final state. It in- 
cludes the angles of the decay products and their masses 
in predefined frames. The decay amplitudes are calcu- 
lated for each event and the "unknown" production am- 
plitudes are varied to obtain the best match between the 
predicted and the observed intensity distribution through 
a maximum likelihood fit. 

The ln(likclihood) function, in its final form, is written 
as: 



ln(£) = £]n 
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The first sum is over the number of events in a given 



4 



mass bin in which a fit is done. The argument of the In is 
just the intensity for each event, /(t,-). The experimen- 
tal acceptance, determined by Monte Carlo simulation, 
was incorporated into the PWA as a normalization fac- 
tor on a wave by wave basis. Two sets of normalization 
integrals were calculated, the accepted, and the raw, 
normalization integrals. The accepted normalization 
integrals were defined as: 



"" M a 



1 M a 
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Where M a is the number of accepted Monte Carlo 
events in a given ujtt mass bin. The raw normalization in- 
tegrals are defined the same, but calculated from the raw 
normalization integrals, over the number of raw events in 
a given lott mass bin. 

With r] — M a /M r , the acceptance corrected intensities 
are written as: 
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Numerous fits were performed on the final sample 
of 168 thousand W7r~ events with various sets of al- 
lowed partial waves in the fit, different lott~ mass bin 
widths, two different regions in —t, different starting val- 
ues for the fit parameters, and different ranks (rank=l 
and rank=2). The general features of the fits did not 
change significantly in any case and since the rank=2 fit 
results did not significantly improve the description of 
the data, a rank=l fit proved to be sufficient. 

The set of waves included in the final PWA fit consisted 
of 21 waves with an lo in the final state, as shown in 
Table |, and 16 waves with ai(1260) and 02 (1320) isobars 
in the final state. The choice of the wave set with 01/02 
isobars was based on obtaining an adequate description 
of the angular distribution of the u sideband events p7j . 
The partial wave intensities for waves with J > 4 were 
found to be insignificant and were eliminated from the 
fits at an early stage of the analysis. A non- interfering 
"isotropic" background wave was included in the fit as a 
cumulative effect of all the small waves omitted from the 
fit. 



TABLE I: Set of 21 waves included in the final PWA fit, with 
an lo in the final state. 
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The acceptance-corrected intensities from various con- 
tributions are shown in Figure ^. The three sets of 
points from top to bottom correspond to contributions 
from waves with an lo in the final state, a non-interfering 
isotropic background wave, and waves with either an a\ 
or an 02 in the final state. Each point in the plots is the 
result of an independent PWA fit for the events in a 60 
MeV wide mass bin. As shown in the figure, the total 
intensity of the 01/02 waves was less than 2% of the to 
waves, and for brevity, they are not listed individually. 
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FIG. 4: Acceptance-corrected intensities from various contri- 
butions: Filled circles: waves with an lo in the final state, 
Squares: "isotropic wave, Open circles: waves with either an 
ai or an 02 in the final state. 



Figure || shows the contributions from the different J pc 
states with an lo in the final state. The major contribu- 
tors to the total intensity are for J pc = 1 H (dominated 
by 6i(1235)), and JP C = 3— (dominated by p 3 (1690)). 
We wish to also note that there appears to be a signifi- 
cant enhancement in the J pc = 2^ intensity at ~ 1650 
MeV. A detailed study of this potentially exotic state will 
be the subject of a future publication. 

The individual partial wave contributions for J pc 
1 with different rrfL are shown in Figure ^ and the 
corresponding phase differences between the D and S 
waves of the same m e are shown in Figure |?]. It is clear 
from the intensity plots that b\ production is dominated 
by the natural parity exchanges. As expected, the phase 
difference for the to c = + and 1 + waves are approxi- 
mately constant, even over the extended b\ mass region 
(~ 1.1 — 1.5 GeV) . The behavior of the to € = 1 ~ phase dif- 
ference (not shown) is erratic due to the small intensities 
and the large error bars. Only the positive reflectivity 
waves are used in the measurement of the D / S. 
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FIG. 5: Acceptance-corrected intensities from the results of a 
PWA fit. The list of ui waves included in this fit is shown in 
Table |l| Individual contributions for different J pc are shown. 
Each J pc is the sum of the allowed m e L included in the fit. 



IV. D/S MAGNITUDE AND PHASE 
MEASUREMENT 



From the final set of PWA fits we determine the in- 
dividual J pc = l" 1 m e production amplitudes, where 
77i c = + ,l ± . As mentioned earlier, the — (+) reflec- 
tivity corresponds to an unnatural(natural) parity ex- 
change. The m £ = + and 1 + b\ production mecha- 
nism is most likely through u> exchange. Since the PWA 
results show that b\ production through the unnatural 
parity exchanges is small and the error bars on the cor- 
responding negative reflectivity waves is large, they were 
omitted from the measurement of the D/S. The D/S 
ratios of both the m e = + , I + b\ decay amplitudes were 
set to a complex number, Re 1 ^ , leaving all other partial 
waves free to vary independently. A grid search was then 
performed in R and <fi for which the -ln(likelihood) func- 
tion, as written in Eq. (|^), was minimized. Convergence 
in both R and 6 was reached after few iterations. The 



details of the procedure can be found in |26|. 

The projections of the -ln(likelihood) function at the 
minimum for D/S magnitude and phase are shown in 
Figures || and || respectively, for the set of waves chosen 
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FIG. 6: Acceptance-corrected intensities for the J pc = 
l + ~m e L partial waves. Only the positive reflectivity waves 
are used in the measurement of D/S. Notice the different 
ordinate scale for the weaker negative reflectivity signals. 



in the PWA fit as shown in Table ffl. These results are 
based on a set of wtt - events in a 160 MeV wide mass bin 
around the b\ mass (1.155—1.315 GeV) with the —t in the 
range (0.1 — 1.5 GeV 2 ). The points in each plot were fit 
to a second order polynomial function where the minima 
are found to be \D/S\ = 0.269 ± (0.009) sta t ± (0.01) sys , 
and <t>{D - S) = 0.184 ± (0.042) stat ± (0.07) sys rad 
or 10.54 ± (2.4) sta t ± (3-9) sys deg. The statistical er- 
ror in each measurement corresponds to the change in 
-ln(likclihood) by 0.5 units. The main sources of the sys- 
tematic error considered were the choice of the wave set 
in the PWA fits and the size of the ujtt~ mass bins used in 
the scans. Regarding the choice of the wave sets, two sets 
of waves with reasonable fit results were used. One set 
consisted of the waves listed in Table | and another con- 
sisted of a subset of 28 largest waves from that list. For 
the u)7r~ mass widths, six sets of independent scans were 
performed in 60, 80, 100, 120, 140 and 160 MeV wide 
wtt mass bins around the nominal b\ mass. No signifi- 
cant systematic change in either \D/S\ or <fi(D — S) was 
observed within the statistical errors, and the systematic 
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FIG. 7: Phase differences between the 1 H m e D and S waves 
with positive reflectivity, m e = + (I + ) in the top (bottom) 
plot. The dashed lines are drawn at the mass of the 6i to 
guide the eyes. 



errors quoted are a conservative estimate. 

Various studies were performed in order to determine 
the leakage from other waves into all the positive reflec- 
tivity J pc = l" 1 waves. In addition, leakage from the 
jvc — \-\ ^ 5*- wave to I?- wave for a given m e = 1~ was 
also investigated. In all cases, the leakage due to experi- 
mental resolutions was determined to be negligible. 



V. CONCLUSIONS 

We have made a study of the reaction n~p — > lutt~ , 
u> — > 7r + 7T~7r° at 18 GeV/c. A partial wave analysis of 
168 thousand events, consistent with the wk~p hypoth- 
esis, shows that the ujtt~ data below 1.6 GeV are domi- 
nated by the J pc = l" 1 resonance 6i(1235). The S-w&ve 
and _D-wave amplitudes of the b\ —> ujit~ decay have been 
determined. It is found that the ratio of the amplitudes, 
\D/S\ = 0.269 ± (0.009) sta t ± (0.01) sys . This represents 
nearly a factor three improvement in error over the cur- 
rent PDG2000 average value of 0.29 ± 0.04. We have 
also determined the phase difference, (f>(D — S) = 0.184± 
(0.042) stat ±(0.07) sys radians, or 10.54±(2.4) stat ±(3.9) sys 
deg, for which no prior results exist. Figure ^| shows that 
the minimum value of <p(D — S) allowed by the above 
errors differs from zero by more than 15a. The level 
of significance varies in other acceptable fits, but it is 
greater than 5a in all cases. 

We can compare our results with the theoretical pre- 
dictions which exist in the literature. Ackleh, Barnes 
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FIG. 8: The projection of the -ln(likelihood) as a function of 
\D/S\. The distribution was fitted to a second order polyno- 
mial function, with the minimum at \D/S\ = 0.269 ± 0.009. 
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FIG. 9: The projection of the -ln(likelihood) as a function 
of 4>{D — S). The distribution was fitted to a second order 
polynomial function, with the minimum at (j>(D—S) = 0.184± 
0.042 rad, or 10.54° ± 2.4°. 



and Swanson Q have made numerical predictions for 
\D/S\ ratios in their 3 Po hadronic decay model as a 
function of the oscillator parameter (3. For the then 
current values of \D / S\ for b\ — > lutt and a\ —> pn, 
they found the best fit value of = 0.448 GeV and 
the corresponding \D/S\ = 0.219. They note, how- 
ever, that these results are in disagreement with "the 
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decay rates of light L(qq) = and L(qq) — 1 mesons 
(which) support — 0.40 GeV". We note that our result, 
\D/S\ = 0.269 ± 0.013 corresponds to /3 = 0.409 ± 0.008, 
in excellent agreement with that for L = and L = 1 
light meson decays. 

Ackleh et al. also make the parameter in- 
dependent prediction that the ratio of ratios, 
R = [{D/S) ai ^ pix ]/[{D/S) bl ^] = -0.5. Us- 
ing the current PDG2000 average, \(D/ S) ai ^ p7I } = 
-0.107 ± 0.016 with our result for (D/S) for b l cjtt, 
we obtain R = —0.40 ± 0.06, with the error determined 
almost entirely by that in the a\ decay ratio. This points 
at the need to improve the measurement of (D/S) for 
ai — ► pit, which is plagued by uncertainty in the Deck 
contribution to this decay. 

The phase difference between the D and S wave de- 
cays of b\ — > ujtt can arise from the final state inter- 



action between the tu and 7r. In a quark interchange 
model calculation, Barnes, Black, and Swanson J?j pre- 
dict <j)(D — S) = 14 deg. Our measurement, <fi(D — S) = 
10.54 ± (2.4) sta t ± (3.9) sys deg, is consistent with their 
prediction. 
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